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ABSTRACT 


The  effect  of  surface  roughness  on  the  radar  backscattering  cross  section 
of  a  perfectly  conducting  nominally  spherical  target  is  examined  by  applying 
the  Kirchhoff  method. 

It  is  shown  that,  for  the  type  of  roughness  and  sphere  size  to  which  the 
Kirchhoff  method  is  applicable,  the  standard  deviation  of  the  cross  section 
increases  with  frequency  according  to  the  law  2V2cr^kt,  until  the  first  Fresnel 
zone  reduces  in  size  to  the  scale  length  of  the  roughness.  At  this  point  a  knee 
in  the  curve  occurs  and  its  further  course  is  determined  by  a  more  detailed 
statistical  description  of  the  surface.  Here  <r  is  the  nominal  cross  section, 

L,  is  the  standard  deviation  of  the  surface  height  h  and  k  =  2tt/\,  where  \  is 
the  wavelength. 


The  average  cross  section  is  shown  to  be  given  by  cr^  1  + 


In 


this  connection,  an  error  that  may  be  significant,  occurring  in  the  work  of 
Hiatt  et  al.on  roughness  effects,  is  pointed  out. 

Some  experimental  results  are  reported  which  support  the  theoretical 
conclusions  and,  moreover,  indicate  that  they  may  be  useful  even  when  the 
scale  length  of  the  roughness  is  smaller  than  the  wavelength. 

Further  theoretical  results  are  included  concerning  the  effect  of  rough¬ 
ness  on  a  general  second  order  surface  and  the  correlation  function  for  the 
cross  section. 


Accepted  for  the  Air  Force 
Stanley  J.  Wisniewski 
Lt  Colonel,  USAF 
Chief,  Lincoln  Laboratory  Office 
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Radar  Echo  Variations  of  a  Large  Rough  Sphere 


I.  INTRODUCTION 

Variation  with  aspect  of  the  radar  backscattering  cross  section  of  a  nom¬ 
inally  spherical  metal  body  detracts  from  its  usefulness  as  a  standard  of 
cross  section.  The  situation  is  deteriorated  further  if  the  mean  cross  section 
is  appreciably  different  from  the  nominal  cross  section.  For  manufacturing 
purposes,  therefore,  it  is  important  to  know  how  the  cross  section  depends 
upon  the  surface  tolerance. 

The  effect  of  roughness  having  a  scale  length  short  compared  with  a  wave¬ 
length  has  been  discussed  by  Hiatt^et  al  (1960)  by  a  method  based  on  an 
impedance  boundary  condition.  Their  theoretical  work  is  limited  to  finding  an 
expression  for  the  average  backscattered  field  strength,  although  useful  ex¬ 
perimental  results  are  presented  concerning  the  variation  of  the  cross  section. 

At  the  other  extreme  of  the  scale  length  of  the  surface  roughness,  that  is, 
for  perturbations  which  smoothly  deform  the  sphere  into  some  still-convex 
body  like  an  ellipsoid,  for  example,  the  problem  is  easily  tackled  by  applying 
the  well  known  relationship  between  the  cross  section  and  the  radii  of  curvature 
at  the  specular  point. 

This  paper  presents  the  results  of  an  application  of  the  Kirchhoff  method. 
The  basic  assumptions  made  are  that  the  perpendicular  distance  from  the  mean 
surface  to  the  true  surface  is  small  compared  with  a  wavelength  and  that  the 
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scale  length  of  the  roughness  is  small  compared  with  the  radius  of  the  sphere 
but  large  compared  with  the  wavelength.  It  is  further  assumed  that  the  sur¬ 
face  is  statistically  uniform  and  isotropic. 

The  scale  length  assumption,  as  it  stands,  limits  the  application  of  the 
results  for  practical  purposes  to  spheres  which  have  had  the  small  scale 
roughness  removed  by  polishing.  However,  experience  with  the  Kirchhoff 
method  has  shown  that  it  yields  accurate  results  even  when  the  scale  length 
assumption  is  relaxed  considerably  (Bechmann  and  Spizzichino,  1963:1V).  In 
Section  VII  some  experimental  results  for  rough  spheres  are  reproduced 
which  tend  to  support  this  view. 

The  Kirchhoff  approach  reduces  the  problem  to  a  form  which  has  already 
received  considerable  attention  in  the  random  phase-screen  theory  of  radio¬ 
star  scintillations  induced  by  the  ionosphere,  [e.g.  Booker,  et  al  (1950), 
Mercier  (1962)]. 

II.  THE  KIRCHHOFF  INTEGRAL 

The  scattering  properties  of  a  perfectly  conducting  sphere  can  be  charac¬ 
terized  by  the  parameter  ka,  where  a  is  the  radius  of  the  sphere  and  k  =  Zu\9 
\  being  the  wavelength  of  the  incident  radiation.  For  large  values  of  ka  the 
backscattered  field  can  be  expressed  as  the  sum  of  an  optics  terms  and  a 
creeping-wave  term.  For  values  of  ka  larger  than  about  10,  the  creeping- 
wave  contribution  is  less  than  one  tenth  that  of  the  optics  term  (Senior,  1964). 
Since,  in  addition,  the  creeping  waves  traverse  a  path  on  the  surface  which 
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is  at  least  ira.  in  length,  it  is  to  be  expected  that  the  effect  of  the  surface  per¬ 
turbations  will  average  out,  in  some  sense,  thereby  making  the  creeping-wave 


term  significantly  less  sensitive  to  the  roughness  than  is  the  optics  term.  It 
is  reasonable  to  suppose,  therefore,  that  the  effect  of  neglecting  the  creeping- 
wave  term  will  not  be  serious  for  even  smaller  values  of  ka  than  10. 

The  optics  term,  to  the  first  order  in  ka,  is  that  obtained  by  the  Kirchhoff 
method  of  calculating  the  scattered  field.  (If  the  spurious  shadow  region  con¬ 
tribution  is  discarded,  the  Kirchhoff  solution  is  accurate  to  at  least  the  third 
order  in  ka.  ) 

When  a  linearly  polarized  plane  wave  of  unit  intensity  propagating  in  the 
negative  z  direction  is  incident  upon  the  rough  sphere,  the  Kirchhoff  solution 
for  the  backscattered  electric  field  at  a  point  in  the  far  zone  can  be  written 
as 


(1) 


A 


Here  R  is  the  distance  from  the  sphere  to  the  field  point,  z(x,  y)  is  the  height 


of  the  illuminated  hemisphere  from  the  plane  z  =  0  and  the  integral  is  carried 


out  over  the  projected  area  of  the  illuminated  hemisphere  on  the  plane  z  =  0. 
A  superfluous  constant  phase  factor  has  been  omitted  and  E  is  understood  to 
be  polarized  in  the  same  direction  as  the  incident  electric  field. 

III.  THE  VARIANCE  OF  THE  CROSS  SECTION 

The  effect  of  the  surface  roughness  may  be  introduced  by  writing  z(x,  y) 
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as 


z(x,  y)  =  z  (x,  y)  +  6(x,  y)  (2) 

o 

where  6  is  the  height  of  the  true  surface,  in  the  z-direction,  from  the  mean 

surface  z  (x,  y).  Thus,  by  definition, 
o 

<6(x,  y)>  =  0  (3) 


where  the  angle  brackets  denote  the  operation  of  taking  the  average  or  expected 
value  of  the  enclosed  quantity,  and  z^(x,  y)  can  be  written 


z  (x,  y)  =V7 

o 


2  2 
x  —  y  —  a 


(4) 


From  (1)  and  (2),  the  expected  value  of  E  is  given  by 


<E>  = 


k 

2ttR 


-  i2kz  (x,  y) 
e  <e 


-  i2k6(x,  y) 


>  dxdy 


(5) 


Figure  1  shows  that,  to  a  first  approximation,  6(x,  y)  is  related  to  h(x,  y)  (the 

height,  in  the  radial  direction,  of  the  true  surface  from  the  mean  surface)  by 

the  equation  6(x,  y)  =  h(x,  y)/cos0,  where  0  is  the  angle  between  the  z  axis  and 

the  radial  direction.  Therefore,  since  the  roughness  is  statistically  uniform 

over  the  surface,  then  in  the  vicinity  of  the  specular  point  <exp{ -i2k6(x,  y)}> 

2  2 

is  a  weak  function  of  the  single  argument  x  +  y  .  In  addition,  a  basic  assump¬ 
tion  here  is  that  h  «  \,  so  that  2k6(x,  y)  is  uniformly  bounded  for  all  k.  This 
clears  the  way  for  a  straightforward  stationary  phase  evaluation  of  (5),  using 


4 


the  form  for  given  by  (4)  with  ka  as  the  large  parameter. 

The  result  is  <E>  =  ia  <exp{—  i2k6(0,  0)}>/2R,  which,  since  6(0,  0)  =  h(0,  0), 
h/\<<  1  and  <h(0,  0)>  =  0,  may  be  written 


<E>  =  (1  -  2k2^2)  i-^- 


(6) 


2  2 

where  =  <h  (x,  y)>,  or,  alternatively. 


<E>  =  (1  -  2k2;2)E 

O 


(7) 


where  E  is  the  backscattered  field  of  the  perfect  sphere  of  radius  a. 
o 

It  is  shown  in  the  Appendix  that,  for  small  perturbations  of  the  scattered 
field,  the  normalized  variance  D{cr}  of  the  cross  section  a  is  given  by 


D{«r)  =  < 


q  —  <q> 
<q> 


>  =  2  < 


E  -  <E> 
<E> 


+  Re 


ffE  -  <E>"12'\ 

}>  * 


(8) 


where  Re  denotes  the  operation  of  taking  the  real  part.  Therefore,  from 
(1),  (2),  (5)and(6),  this  may  be  written,  to0{(k£)^},  as 


D{q}  = 


2k 


2  2 
it  a 


[I 


-i2k[  z  (x,  y)  -  z  (x\  y1)] 


F^(x,  y,  x'f  y1)  dxdydx'dy' 


—  Re  <  \  e 


-i2k[  zQ(x»  y)  +  Zq(x',  y ’ ) ] 


F^(x,  y,  x',  y')  dxdy  dx'dy,'j- 


.  (9) 
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where 


F,  u,  y,  x\  y')  =  <e’“k[s(x*  y>-  6<x''  y'>]>_<e-i2k6<*.  y)><ei^(x'.  y')> 


F,  (X.  y.  x',  y')  =  Ce’1*1  6(x>  y)  +  5(x’’  y')]>  -  <e -i2k6|x’  y)>  <*-»»«(*'.  v\ 


It  is  clear  that  both  and  F^  are  zero  when  the  distance  between  the 
points  (x,  y)  and(x',  y 1 )  is  large  enough  for  the  corresponding  6!s  to  be  uncor¬ 
related.  Therefore,  when  the  scale  length  L  of  the  surface  roughness  is 
smaller  than  the  width  w  of  the  first  Fresnel  zone,  the  variation  of  F^  and  F^ 
as  functions  of  x  —  x',  y  —  y1  is  stronger  than  that  of  the  exponential  factors 
and  the  method  of  stationary  phase  is  not  directly  applicable. 

However,  the  period  of  the  exponential  factors  in  the  integrands  in  (9) 
steadily  increases  as  x,  x’,  y  or  y1  increase,  so  that  eventually  the  period  is 
small  compared  with  the  scale  length  of  the  roughness.  Such  regions  contrib¬ 
ute  negligibly  to  the  integrals,  which  leads  to  the  conclusion  that  the  effective 
region  of  integration  in  x,  x’,  y,  y1  space  is  a  hyper-sphere  centered  at  the 


— 2 - 2 

x  +  y  at  which  the 


o  o 

period  of  the  exponential  factors  is  equal  to  the  scale  length  L 


Vi--2.  2-.  2 


(x  +  y  )/a 

of  the  roughness.  (The  scale  length  on  the  x,  y  plane  is  a  projection  of  the 

true  surface  scale  length.  )  The  value  of  r  is  easily  deduced  to  be  a\/4L». 

o 

Since,  by  assumption,  X/L  «  1,  r^  is  small  compared  with  the  radius  of  the 
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sphere. 


Therefore,  for  all  L,  the  significant  contribution  to  the  integrals  in  (9) 

derives  from  a  small  region  about  the  origin  within  which  6(x,  y)  and  z  (x,  y) 

2  2 

are  given  essentially  by  h(x,  y)  and-(x  +  y  )/2a,  respectively.  In  addition, 
since  the  limits  of  integration  do  not  lie  within  or  close  to  this  region,  they 
may  be  extended  to  infinity.  D{<r}  can  then  be  written  as 


.  k  r  2  2  ,2  ,2, 

1  7  lx  +  y  -  x'  -  y1  J 


(x  —  x1)2  +  (y  —  y1)2  J  dxdx1  dydy1 


■{f 


+  Rei  \  e 


•  k  r  2  4.  2  ^  ,2  4.  ,2i 

l  —  [  x  +  y  +  x‘  +  y  J 

a 


(x  —  x1  )2  +  (y  ~  y’)2)  dxdx1  dydy'j* 


,(10) 


where  ^2Cu(~/(x  —  x')2  +  (y  —  y')2  )  =  <h(x,  y)  h(x',  y')>,  (3^(0)  =  1  and  the  fact 

2  2 

that  h/\  <<  1  has  been  used  to  reduce  F  and  F  to  the  expressions  4k  g  C 

12  h 

and  —4k  £,  C  respectively.  As  defined  here  C  is  the  correlation  function 
h  h 


for  the  surface. 

By  means  of  the  substitutions  x  =  u  +  v/ 2,  x1  =  u  —  v/2,  y  =  u1  +  v1  /2, 


yi  -  ui  _  v' / 2,  the  form  for  D{cr}  given  by  (10)  simplifies  to 
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UJ 

D{o-}  =  8k\  ^  1  —  ^  C^C,2ap/k)  sinp  dp 


(ID 


This  is  the  required  result. 

When  the  scale  length  is  small  compared  with  the  width  w  of  the  first 


Fresnel  zone  (w  =  V 27Ta/k),  the  integral  in  (11)  is  clearly  small  compared 
with  unity.  Thus  in  this,  the  case  of  most  practical  interest,  the  normalized 
standard  deviation  of  the  cross  section  is  given  simply  by  2V2 k£. 

IV.  THE  AVERAGE  CROSS  SECTION 

2 

The  cross  section  is  proportional  to  I  E  1  ,  so  that  the  average  normalized 
cross  section  <(r>/or^  is  given  by 


<cr> 

0" 


<|E|2> 


I  <E>  I 


(12) 


The  rearrangement  of  the  right  side  includes  the  term  <|(E  —  <E>)/<E>|  > 

which  has  been  evaluated  in  Section  III.  It  is,  from  (8),  just  half  the  first 

term  on  the  right  of  (11),  or  4k^C,^.  Therefore,  (12)  can  be  rewritten,  using 

2  2  ^  2  2 

(7),  as  <c r>/cr  =  (1  —  2k  t,  )  (1  +  4k  £,  ),  which  is  unity  to  the  order  in  kh  to 

o 

which  the  approximation  has  been  carried.  That  is 


1+  O-j(kh)3} 


(13) 
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This  result  leads  one  to  conclude  that  a  departure  of  the  average  cross 
section  from  the  nominal  cross  section  is  unlikely  to  be  measureable  under 
the  conditions  for  which  the  theory  is  valid. 

It  is  of  interest  at  this  point  to  mention  the  work  of  Hiatt  et  al  (I960),  who 
use  an  impedance  boundary  condition  to  represent  small  scale  surface  rough¬ 
ness.  They  find  an  expression  for  the  expected  value  of  the  backscattered 

field  and  then  square  this  to  obtain  the  expected  value  of  the  backscattering 

Z  Z 

cross  section.  This  is  clearly  in  error,  for  <  |  E  |  >  =  <  | <E>  +  (E  —  <E>)  |  > 

2  2 

=  |  <E>  |  +  <IE  —  <E>  I  >,  so  that  they  have  neglected  to  include  the  term 

2 

<  |  E  —  <E>  |  >. 

The  application  of  the  Kirchhoff  method  to  the  same  problem  has  shown 
that  when  this  additional  term  is  included  it  just  cancels  the  effect  of  the 
roughness  on  the  term  |<E>|^.  Hence,  although  the  two  approaches  to  the 
problem  do  not  start  from  the  same  assumptions,  the  Kirchhoff  method  indi¬ 
cates  the  necessity  of  including  the  additional  term.  One  concludes,  therefore, 
that  a  modification  of  the  average  backscattering  cross  section  due  to  small 
scale  roughness  has  yet  to  be  demonstrated  theoretically. 

V.  GENERALIZATIONS 

It  is  a  simple  matter  to  rework  the  derivation  for  a  mean  surface  in  the 
form  of  a  general  second-order  surface.  Instead  of  the  single  radius  a,  one 
has  to  deal  with  the  separate  principle  radii  of  curvature  a^  and  a^  at  the 
specular  point. 
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In  place  of  (11),  one  finds 


2tt  oo 


D{cr}  =  8kZ;Z 


1  —  — ~  ^  C^V 2 (a ~ cos ^0  +  a^sin^Jp  /kl  sinp  dp  d0 


o  o 


so  that  again,  when  the  scale  length  is  small  compared  with  the  least  width  of 
the  first  Fresnel  zone,  the  integral  is  small  and  the  normalized  standard  de¬ 
viation  of  the  cross  section  is  2V2k£>. 

Equation  (13)  for  the  average  cross  section  is  found  to  apply  unaltered  to 
the  more  general  surface. 

A  further  result  of  some  interest  is  the  correlation  function  C  (0)  for  the 

<T 

backscattering  cross  section,  which  may  be  defined  as 


C  (0) 
0" 


R(9) 

R(0) 


(14) 


where  R(0)  =  <(r(0  )  o-(0^+  0)>  —  <cr> ^  and  the  arguments  0  and  0^+  0  specify 
the  directions,  based  on  an  arbitrary  great  circle  on  the  surface  of  the  sphere, 
at  which  the  cross  sections  are  measured. 

The  derivation  is  straightforward  provided  a0  is  small  compared  with  the 
radius  of  the  sphere,  for  then  it  is  sufficient  to  replace  x’  and  y'  in  the  ex¬ 
ponential  factors  in  the  integrands  in  (9)  by  x'  —  x^  and  y’  —  y  ,  where 
2  2  2  2 

x  +  y  =  a  0  .  The  result  may  be  expressed  as 
o  o 


10 


CO 

Jo(0^2kap)sin(^p  +  ka0^/2^dp  , 

o 

(15) 

where  J  (u)  is  the  Bessel  function  of  the  first  kind  of  order  n  and  argument  u. 
n 

When  the  scale  length  of  the  surface  perturbations  is  small  compared 

with  the  width  of  the  first  Fresnel  zone,  the  integral  in  (15)  is  clearly  small 

compared  with  the  first  term  in  brackets.  In  this  case  (15)  reduces  to 
2  2  2 

R(0)  =  8k  t  cr  ~C  (a0),  so  that  the  correlation  function  for  the  cross  section 
o  h 

is,  from  (14),  given  simply  by  C  (0)  =  C^(a0). 

VI,  NUMERICAL  STUDY 

The  derivation  of  expression  (11)  for  the  variance  of  the  cross  section  is 
based  on  the  well  tried  Kirchhoff  method.  Thus  if  the  result  is  incorrect  it  is 
likely  to  be  because  the  derivation  is  faulty  rather  than  the  Kirchhoff  method. 

To  provide  a  check  on  the  derivation,  a  computer  program  was  written 
to  compute,  using  the  Kirchhoff  method,  the  backscattering  cross  section  for 
a  sphere  with  a  particular  type  of  random  surface.  By  repeating  the  compu¬ 
tation  for  a  number  of  different  samples  of  the  same  random  process  the 
variance  of  the  cross  section  was  obtained. 

The  random  surface  was  generated  by  dividing  the  surface  of  the  mean 
sphere  into  a  large  number  of  curvilinear  squares  of  roughly  equal  area.  The 
true  surface  was  then  defined  to  be  given  in  spherical  coordinates  by 


R(0 )  =  8k Z;2o-  Z 
o 


C, (a0) 
h 


"I  ch(^ 


/k 
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r(0,  <J>)  =  a  +  h(0,  4>),  with  h(0,  <j>)  being  constant  for  all  points  (0,  <j))  lying 
within  a  particular  curvilinear  square.  The  value  of  h  assigned  to  each  square 
was  an  independent  sample  of  a  Gaussian  random  process  of  zero  mean  and 
given  variance. 

The  results  are  shown  in  Table  I.  The  entries  in  columns  1  —  5  give  the 
various  ratios  of  significance  to  the  computation,  with  L  being  taken  as  twice 
the  side  length  of  the  elemental  surface  areas.  The  entries  in  column  6  are 
the  corresponding  values  of  the  normalized  standard  deviation  of  the  cross 
section  given  by  the  expression  Z\l  2k£,,  which,  according  to  the  derivation  is 
appropriate  for  small  values  of  L/w.  The  normalized  standard  deviation  of 
the  N  cross  section  values  obtained  with  the  computer  program  is  entered  in 
column  7  and  N  is  recorded  in  column  8. 

TABLE  I 

COMPUTATION  OF  THE  NORMALIZED  STANDARD 
DEVIATION  OF  CROSS  SECTION  FOR  A  ROUGH  SPHERE 


1 

2 

3 

4 

5 

6 

7 

8 

ka 

4A 

L/\ 

L/a 

(L/w)2 

2'v/2kt, 

\/D{o-} 

N 

50 

0.  8  x  10-3 

2.  2 

0.  28 

0.  65 

0.  014 

0.  017 

19 

100 

1.6x10  3 

4.  6 

0.  28 

1.3 

0.  028 

0.  029 

19 

200 

3.  2  x  10"3 

9.  1 

0.  28 

2.  6 

0.  056 

0.  031 

19 

200 

3.  2  x  10"3 

2.  2 

0.  072 

0.  16 

0.  056 

0.  051 

12 

12 


The  first  three  rows  show  the  effect  of  increasing  the  frequency,  leaving 
all  else  unchanged.  The  agreement  between  columns  6  and  7  is  good  except 
at  the  highest  frequency,  where,  in  view  of  the  fact  that  (L/w)^  has  the  value 
2.6,  it  is  to  be  expected  that  the  simple  form  Z^lZkL,  should  over-estimate 
*vD {cr}’.  The  final  row  shows  the  effect  of  then  reducing  L  by  a  factor  of  0.  25, 
which  makes  (L/w)^  small  again  and  the  agreement  is  regained. 

VII.  EXPERIMENTAL  RESULTS 

The  theoretical  work  reported  here  was  carried  out  as  a  part  of  the  pro¬ 
gram  which  culminated  in  the  successful  launch  of  a  113  cm  radar  calibration 
sphere  (Object  number  1361)  on  May  6,  1965  into  a  nearly  circular  1500 
nautical-milt  orbit  of  32.  11°  inclination.  The  initial  attempt  to  manufacture 
the  sphere  involved  the  welding  together  of  two  hemispherical  aluminum  shells. 
It  was  found  that  the  welding  operation  distorted  the  shells  over  a  wide  equa¬ 
torial  zone  to  a  degree  that  made  the  resulting  sphere  unacceptable  as  a  radar 
calibration  device.  For  this  reason,  a  different  method  of  joining  the  two 
hemispheres  was  used  for  the  flight  sphere.  However,  the  welded  sphere  was 
useful  for  surface  roughness  experiments  because  the  variations  in  cross 
section  were  large  enough  to  enable  good  measurements  to  be  made  on  a 
standard  radar  cross-section  range. 

The  sphere  was  polished  after  welding,  so  that  the  final  condition  of  the 
surface  was  very  smooth  with  undulations  of  the  order  of  0.  1  cm  in  height 
and  15  cm  in  length.  The  surface  was  not  statistically  uniform,  since  the 
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surface  was  most  "rough"  in  the  region  of  the  equator.  Therefore,  since 

this  immediately  violates  one  of  the  conditions  necessary  for  the  simple  theory, 

it  was  decided  that  a  detailed  dimensional  survey  of  the  surface  would  not  be 

worthwhile.  Accordingly,  the  two  statistical  characterizations,  £,  and  C  (p), 

h 

of  the  surface  necessary  to  evaluate  Dio-}  using  (11)  were  estimated  from  a 
single  sectional  profile  of  the  sphere  taken  by  means  of  an  accurate  circular 
check-ring  lying  in  the  equatorial  plane  (the  plane  of  the  weld  line).  The  dis¬ 
tance  from  the  check-ring  to  the  surface  was  measured  at  10°  increments 
round  the  full  circle. 

The  check-ring  measurements  were  reduced  to  equivalent  "displacement- 

from-mean-sphere"  measurements  by  subtracting  from  them  their  mean  and 

their  first  Fourier  component.  The  standard  deviation  £,  of  the  reduced 

measurements  was  found  to  be  0.  0591  cm  and  the  correlation  function  (defined 

to  be  unity  at  0°)  was  found  to  be  0.  607  at  10°  and  -0.  017  at  20°.  The  last 

number  indicates  that  the  correlation  function  falls  quickly  towards  zero  and 

in  addition,  since  the  surface  is  smooth,  the  slope  of  the  correlation  function 

must  be  zero  at  the  origin.  This  indicates  that  an  appropriate  simple  model 

2  2 

for  the  correlation  function  is  C^(p)  =  exp(-p  /L  )  with  L  =  13.  95,  the  value 

required  to  make  C  (p)  equal  to  0.  607  when  p  is  the  surface  arc  length  equiv- 

h 

alent  to  an  angular  displacement  of  10°. 

Substituting  this  form  for  C^(p)  into  (11)  leads  to  the  following  expression 
for  the  standard  deviation  of  the  backscattering  cross  section. 


14 


Vdm  = 


2 


(16) 


Z\jzk^ 


7T- 


W 


» 


where  w  =  \Za\  is  the  width  of  the  first  Fresnel  zone. 

Figure  2  is  a  graph  of  expression  (14)  as  a  function  of  frequency  using  the 
values  of  £,  and  L  obtained  as  described  from  the  check-ring  measurements 
and  the  nominal  sphere  radius  of  56.  5  cm.  The  broken  straight  line  is  the 
low  frequency  asymptote. 

The  measured  values  of  the  standard  deviation  of  the  backscattering  cross 
section  are  entered  at  points.  Measurements  were  made  at  UHF  and  in  the 
bands  L,  S,  X,  and  K.  The  plotted  points  were  computed  from  the  360°  range 
record  taken  around  the  equator  or  weld  line  of  the  sphere. 

In  Table  II  the  theoretical  and  experimental  data  are  presented 
numerically. 
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TABLE  II 


COMPARISON  OF  THEORY  AND 


EXPERIMENT  FOR 


F  requency 
(GHz) 

ka 

W  avelength 
(cm) 

0.  4 

4.  73 

75 

1.  2 

14.  2 

25 

3.  3 

39.  1 

9.09 

8.  5 

100.  6 

3.  53 

24.  0 

284 

1.  25 

113  cm  SPHERE 
\[D{cr} 

Theoretical  Measured  Measured  (db) 


0. 0137 

0.018 

0.  078 

0. 0385 

0.  046 

0.  20 

0. 0742 

0.  064 

0.  28 

0. 0923 

0.  087 

0.  38 

0.  0965 

0.  187 

0.  86 

The  standard  deviation  of  the  measured  values  is  also  expressed  in  decibels 
for  convenience. 

The  agreement  between  theory  and  experiment  is  surprisingly  good  in 
view  of  the  approximations  used  to  describe  the  surface.  The  fact  that  at 
24  GHz,  the  agreement  is  lost  is  probably  due  to  the  coarse  measuring  interval 
(10°)  used  to  plot  the  sectional  profile.  For  at  24  GHz  the  first  Fresnel  zone 
is  about  8.  5°  wide,  so  at  this  frequency,  smaller  scale  roughness,  undisclosed 
by  the  coarse  measuring  interval,  comes  into  play. 

It  should  be  mentioned  that  at  the  three  higher  frequencies  it  was  possible 
to  expand  the  scale  factor  of  the  output  pen  recording  to  a  value  typically  in 
the  region  of  1  db/in.  The  equipment  used  at  0.  4  and  1.  2  GHz,  however,  was 
not  amenable  to  this  modification  with  the  result  that  the  scale  factor 
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remained  at  about  6  db/in.  Thus  the  figure  of  0*  078  db  for  the  standard  devi¬ 
ation  of  the  cross  section  at  0.  4  GHz  was  obtained  from  a  pen  recording  with 
a  peak-to-peak  variation  of  about  0.  05  inch.  In  addition,  the  background 
level  during  the  measurement  was  of  the  order  of  -40  db,  which  is  too  high  to 
be  negligible.  Therefore,  the  apparent  agreement  at  0.  4  GHz  is  likely  to  be 
more  fortuitous  than  real. 

It  is  interesting  to  note  that  at  the  two  lowest  measurement  frequencies 
the  wavelength  was  larger  than  the  scale  length  of  the  roughness.  This  vio¬ 
lates  one  of  the  conditions  required  by  the  Kirchhoff  method,  but  there  is  no 
apparent  divergence  between  theory  and  experiment  as  the  wavelength  in¬ 
creases. 

Some  further  experimental  data  of  interest  are  contained  in  the  paper  by 
Hiatt  et  al.(1960).  They  measured  the  standard  deviation  of  the  cross  section 
of  a  rough  26  cm  aluminum  sphere  at  two  wavelengths,  3.  1  cm  and  1. 3  cm, 
and  obtained  the  figures  of  0.4  db  and  1.04  db,  respectively.  The  standard 
deviation  of  the  true  surface  from  the  mean  surface  was  estimated  as  0.  037  cm 
and  the  scale  length  of  the  roughness  as  0.  101  cm.  Thus  at  both  frequencies 
the  scale  length  to  wavelength  ratio  was  less  than  one  tenth,  which  places  the 
character  of  the  roughness  well  away  from  the  regime  of  applicability  of  the 
Kirchhoff  method.  However,  a  formal  application  of  (11)  to  this  case  yields 
the  figures  of  0.  92  db  and  2.  2  db  for  the  standard  deviation  of  the  cross 
section  at  the  wavelengths  3.  1  cm  and  1.  3  cm,  respectively.  (For  such  a 
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small  scale  length,  the  abbreviated  form  of  (11)  applies,  making  it  unneces¬ 
sary  to  know  the  form  of  the  correlation  function  for  the  roughness.  ) 

A  second  assumption  necessary  to  derive  (11)  was  that  the  maximum 
perturbation  of  the  cross  section  be  small.  This,  too,  is  not  true  for  the 
present  case,  and  additional  error  is  incurred  in  assuming  a  linear  relation 
between  the  standard  deviation  expressed  as  a  pure  number  and  the  same 
quantity  when  expressed  in  decibels.  In  spite  of  this,  the  theory  agrees  with 
experiment  to  within  a  factor  of  about  two  at  both  measurement  frequencies. 
This  suggests  that  (11)  may  be  a  useful  guide  to  the  effect  of  roughness  even 
for  scale  lengths  considerably  smaller  than  a  wavelength,  but  further  mea¬ 
surements  using  smaller  amplitude  roughness  are  necessary  to  find  out 
whether  this  is  so. 

There  was  no  possibility  of  measuring  the  average  departure  of  the  cross 
section  from  nominal,  because  the  sphere  itself  was  at  least  as  good,  for  the 
purposes  of  calibration,  as  any  other  available  standard. 

VIII.  CONCLUSIONS 

The  effect  of  surface  roughness  on  the  radar  backscattering  cross  section 
of  a  large  metal  sphere  has  been  examined  by  applying  the  Kirchhoff  method. 
The  roughness  has  been  assumed  to  be  statistically  uniform  and  isotropic,  and 
the  effect  on  the  cross  section  has  been  determined  by  obtaining  expressions 
for  the  average  cross  section  and  the  standard  deviation  of  the  cross  section. 

It  is  difficult  to  prepare  spheres  having  rough  surfaces  with  accurately 
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known  and  uniform  statistics,  and  equally  difficult  to  carry  out  accurate  mea¬ 


surements  on  the  small  perturbations  of  the  cross  section,  but  those 
experimental  results  which  have  been  obtained  appear  to  confirm  the  validity 
of  the  theory  within  its  region  of  applicability  and  to  suggest  that  the  theory 
may  be  a  useful  guide  even  for  roughness  of  scale  size  small  compared  a 
wavelength. 
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APPENDIX 


This  appendix  records  the  derivation  of  (8)  in  the  main  text. 

The  backscatter ed  complex  field  E  may  be  written  as 

E  =  <E>  +  6E  (A-l) 

and  if  it  be  assumed  that  |  6E  |  <<  |E|,  then  the  absolute  value  JEJ  of  E  is 
given  closely  by  the  absolute  value  |<E>|  of  <E>  together  with  the  component 
of  6E  in  the  direction  of  <E>.  That  is 


(A-  2) 


where  the  star  denotes  complex  conjugate. 

From  (A-l)  it  is  clear  that  <6E>  =  0,  so  that  by  taking  the  expected  value 
of  (A-2)  one  finds  <|E|>  «  |<E>|  which  enables  one  to  rewrite  (A-2)  as 


Therefore, 


lEj  -  <  I  E  '  > 

<  I  E  !  > 


(A-3) 


since 


<E>  <E>*  and  6E  =  E  —  <E>. 
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Now  the  cross  section,  cr,  is  proportioned  to  1 E I  f  so  that 


D{cr}  =  < 


(T  —  <(T> 
<(T> 


>  =  < 


r  |ej_ 

L  < 


b-  <Ie,2>'1 


(A-4) 


The  notation  6|Ei=  !E|  —  <|E|>  allows  the  right  side  of  (A-4)  to  be  ex¬ 
panded  and  written  as 


IeI2-<IeIV2 


<|E|2> 


>  =  4 


<6 'El  > 
<  I  E  I  >2 


1  +  0 


f  6  I E  1 

t  !EU 


(A-5) 


so  that,  since  both  6|E!  <  |  6E  |  and,  by  assumption,  |  6E  |  <<  )  E  I ,  the  factor 
in  square  brackets  on  the  right  of  (A-5)  is  essentially  unity  and 


D{<r}  «  4 


<6  I  E  I  > 


<  I  E  |  > 


4< 


2 

> 


(A-6) 


By  combining  (A-3)  with  (A-6),  the  required  expression  (8)  is  obtained. 
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/ 


Fig.  1.  Geometry  of  rough  sphere. 


FREQUENCY  (GHz) 


Fig.  2.  Standard  deviation  of  echo  area  as  a  function  of 
frequency  -  theory  and  measurement. 
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